Exploring hypotheses of the actions of TGF-beta 1 in epidermal wound healing using a 3D computational multiscale model of the human epidermis

In vivo and in vitro studies give a paradoxical picture of the actions of the key regulatory factor TGF-beta 1 in epidermal wound healing with it stimulating migration of keratinocytes but also inhibiting their proliferation. To try to reconcile these into an easily visualized 3D model of wound healing amenable for experimentation by cell biologists, a multiscale model of the formation of a 3D skin epithelium was established with TGF-beta 1 literature-derived rule sets and equations embedded within it. At the cellular level, an agent-based bottom-up model that focuses on individual interacting units ( keratinocytes) was used. This was based on literature-derived rules governing keratinocyte behavior and keratinocyte/ECM interactions. The selection of these rule sets is described in detail in this paper. The agent-based model was then linked with a subcellular model of TGF-beta 1 production and its action on keratinocytes simulated with a complex pathway simulator. This multiscale model can be run at a cellular level only or at a combined cellular/subcellular level. It was then initially challenged ( by wounding) to investigate the behavior of keratinocytes in wound healing at the cellular level. To investigate the possible actions of TGF-beta 1, several hypotheses were then explored by deliberately manipulating some of these rule sets at subcellular levels. This exercise readily eliminated some hypotheses and identified a sequence of spatial-temporal actions of TGF-beta 1 for normal successful wound healing in an easy-to-follow 3D model. We suggest this multiscale model offers a valuable, easy-to-visualize aid to our understanding of the actions of this key regulator in wound healing, and provides a model that can now be used to explore pathologies of wound healing

Delayed Cutaneous Wound Healing and Aberrant Expression of Hair Follicle Stem Cell Markers in Mice Selectively Lacking Ctip2 in Epidermis

This is the publisher’s final pdf. The published article is copyrighted by PLoS and can be found at: http://www.plosone.org/home.action.Background: COUP-TF interacting protein 2 [(Ctip2), also known as Bcl11b] is an important regulator of skin homeostasis, and is overexpressed in head and neck cancer. Ctip2(ep-/-) mice, selectively ablated for Ctip2 in epidermal keratinocytes, exhibited impaired terminal differentiation and delayed epidermal permeability barrier (EPB) establishment during development, similar to what was observed in Ctip2 null (Ctip2(-/-)) mice. Considering that as an important role of Ctip2, and the fact that molecular networks which underlie cancer progression partially overlap with those responsible for tissue remodeling, we sought to determine the role of Ctip2 during cutaneous wound healing. \ud \ud Methodology/Principal Findings: Full thickness excisional wound healing experiments were performed on Ctip2(L2/L2) and Ctip2(ep-/-) animals per time point and used for harvesting samples for histology, immunohistochemistry (IHC) and immunoblotting. Results demonstrated inherent defects in proliferation and migration of Ctip2 lacking keratinocytes during re-epithelialization. Mutant mice exhibited reduced epidermal proliferation, delayed keratinocyte activation, altered cell-cell adhesion and impaired ECM development. Post wounding, Ctip2(ep-/-) mice wounds displayed lack of E-Cadherin suppression in the migratory tongue, insufficient expression of alpha smooth muscle actin (alpha SMA) in the dermis, and robust induction of K8. Importantly, dysregulated expression of several hair follicle (HF) stem cell markers such as K15, NFATc1, CD133, CD34 and Lrig1 was observed in mutant skin during wound repair. \ud \ud Conclusions/Significance: Results confirm a cell autonomous role of keratinocytic Ctip2 to modulate cell migration, proliferation and/or differentiation, and to maintain HF stem cells during cutaneous wounding. Furthermore, Ctip2 in a non-cell autonomous manner regulated granulation tissue formation and tissue contraction during wound closure

Development of a Three Dimensional Multiscale Computational Model of the Human Epidermis

Transforming Growth Factor (TGF-β1) is a member of the TGF-beta superfamily ligand-receptor network. and plays a crucial role in tissue regeneration. The extensive in vitro and in vivo experimental literature describing its actions nevertheless describe an apparent paradox in that during re-epithelialisation it acts as proliferation inhibitor for keratinocytes. The majority of biological models focus on certain aspects of TGF-β1 behaviour and no one model provides a comprehensive story of this regulatory factor's action. Accordingly our aim was to develop a computational model to act as a complementary approach to improve our understanding of TGF-β1. In our previous study, an agent-based model of keratinocyte colony formation in 2D culture was developed. In this study this model was extensively developed into a three dimensional multiscale model of the human epidermis which is comprised of three interacting and integrated layers: (1) an agent-based model which captures the biological rules governing the cells in the human epidermis at the cellular level and includes the rules for injury induced emergent behaviours, (2) a COmplex PAthway SImulator (COPASI) model which simulates the expression and signalling of TGF-β1 at the sub-cellular level and (3) a mechanical layer embodied by a numerical physical solver responsible for resolving the forces exerted between cells at the multi-cellular level. The integrated model was initially validated by using it to grow a piece of virtual epidermis in 3D and comparing the in virtuo simulations of keratinocyte behaviour and of TGF-β1 signalling with the extensive research literature describing this key regulatory protein. This research reinforces the idea that computational modelling can be an effective additional tool to aid our understanding of complex systems. In the accompanying paper the model is used to explore hypotheses of the functions of TGF-β1 at the cellular and subcellular level on different keratinocyte populations during epidermal wound healing

Identification of Pax6-Dependent Gene Regulatory Networks in the Mouse Lens

Lineage-specific DNA-binding transcription factors regulate development by activating and repressing particular set of genes required for the acquisition of a specific cell type. Pax6 is a paired domain and homeodomain-containing transcription factor essential for development of central nervous, olfactory and visual systems, as well as endocrine pancreas. Haploinsufficiency of Pax6 results in perturbed lens development and homeostasis. Loss-of-function of Pax6 is incompatible with lens lineage formation and results in abnormal telencephalic development. Using DNA microarrays, we have identified 559 genes expressed differentially between 1-day old mouse Pax6 heterozygous and wild type lenses. Of these, 178 (31.8%) were similarly increased and decreased in Pax6 homozygous embryonic telencephalon [Holm PC, Mader MT, Haubst N, Wizenmann A, Sigvardsson M, Götz M (2007) Loss- and gain-of-function analyses reveals targets of Pax6 in the developing mouse telencephalon. Mol Cell Neurosci 34: 99–119]. In contrast, 381 (68.2%) genes were differently regulated between the lens and embryonic telencephalon. Differential expression of nine genes implicated in lens development and homeostasis: Cspg2, Igfbp5, Mab21l2, Nrf2f, Olfm3, Spag5, Spock1, Spon1 and Tgfb2, was confirmed by quantitative RT-PCR, with five of these genes: Cspg2, Mab21l2, Olfm3, Spag5 and Tgfb2, identified as candidate direct Pax6 target genes by quantitative chromatin immunoprecipitation (qChIP). In Mab21l2 and Tgfb2 promoter regions, twelve putative individual Pax6-binding sites were tested by electrophoretic mobility shift assays (EMSAs) with recombinant Pax6 proteins. This led to the identification of two and three sites in the respective Mab21l2 and Tgfb2 promoter regions identified by qChIPs. Collectively, the present studies represent an integrative genome-wide approach to identify downstream networks controlled by Pax6 that control mouse lens and forebrain development

Using Evolutionary Conserved Modules in Gene Networks as a Strategy to Leverage High Throughput Gene Expression Queries

Background: Large-scale gene expression studies have not yielded the expected insight into genetic networks that control complex processes. These anticipated discoveries have been limited not by technology, but by a lack of effective strategies to investigate the data in a manageable and meaningful way. Previous work suggests that using a pre-determined seednetwork of gene relationships to query large-scale expression datasets is an effective way to generate candidate genes for further study and network expansion or enrichment. Based on the evolutionary conservation of gene relationships, we test the hypothesis that a seed network derived from studies of retinal cell determination in the fly, Drosophila melanogaster, will be an effective way to identify novel candidate genes for their role in mouse retinal development. Methodology/Principal Findings: Our results demonstrate that a number of gene relationships regulating retinal cell differentiation in the fly are identifiable as pairwise correlations between genes from developing mouse retina. In addition, we demonstrate that our extracted seed-network of correlated mouse genes is an effective tool for querying datasets and provides a context to generate hypotheses. Our query identified 46 genes correlated with our extracted seed-network members. Approximately 54% of these candidates had been previously linked to the developing brain and 33% had been previously linked to the developing retina. Five of six candidate genes investigated further were validated by experiments examining spatial and temporal protein expression in the developing retina. Conclusions/Significance: We present an effective strategy for pursuing a systems biology approach that utilizes an evolutionary comparative framework between two model organisms, fly and mouse. Future implementation of this strategy will be useful to determine the extent of network conservation, not just gene conservation, between species and will facilitate the use of prior biological knowledge to develop rational systems-based hypotheses

The human keratins: biology and pathology

The keratins are the typical intermediate filament proteins of epithelia, showing an outstanding degree of molecular diversity. Heteropolymeric filaments are formed by pairing of type I and type II molecules. In humans 54 functional keratin genes exist. They are expressed in highly specific patterns related to the epithelial type and stage of cellular differentiation. About half of all keratins—including numerous keratins characterized only recently—are restricted to the various compartments of hair follicles. As part of the epithelial cytoskeleton, keratins are important for the mechanical stability and integrity of epithelial cells and tissues. Moreover, some keratins also have regulatory functions and are involved in intracellular signaling pathways, e.g. protection from stress, wound healing, and apoptosis. Applying the new consensus nomenclature, this article summarizes, for all human keratins, their cell type and tissue distribution and their functional significance in relation to transgenic mouse models and human hereditary keratin diseases. Furthermore, since keratins also exhibit characteristic expression patterns in human tumors, several of them (notably K5, K7, K8/K18, K19, and K20) have great importance in immunohistochemical tumor diagnosis of carcinomas, in particular of unclear metastases and in precise classification and subtyping. Future research might open further fields of clinical application for this remarkable protein family